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C^ , Context. AUMic is a young M-type star surrounded by an edge-on optically thin debris disk that shares many common observational 
properties with the disk around /JPictoris. In particular, the scattered light surface brightness profile falls off as ~ r^^ outside 120 AU for 
P Pictoris and 35 AU for AU Mic. In both cases, the disk color raises with increasing distance beyond these reference radii. 
Aims. We present the first comprehensive analysis of the AUMic disk properties since the system was discovered by Kalas et al. (2004). We 
explore whether the dynamical model, successful to reproduce the y6 Pic brightness profile (e.g. Augereau et al. 2001), could apply to AUMic. 
/^^ Methods. We calculate the surface density profile of the AU Mic disk by performing the inversion of the near-IR and visible scattered light 
■~>«^ brightness profiles measured by Liu (2004a) and Krist et al. (2005), respectively. We discuss the grain properties by analysing the blue color of 
1-^ the disk in the visible (Krist et al. 2005) and by fitting the disk spectral energy distribution. We finally evaluate the radiation and wind forces 
*T^, on the grains. The impact of the recurrent X-ray and UV-flares on the dust dynamics is also discussed. 

O ' Results. We show that irrespective of the mean scattering asymmetry factor of the grains, most of the emission arises from an asymmetric, 
V^ ' coUisionally-dominated region that peaks close to the surface brightness break around 35 AU. The elementary scatterers at visible wavelengths 
C/j are found to be sub-micronic, but the inferred size distribution underestimates the amount of large grains, resulting in too low sub-millimeter 
Cu emissions compared to the observations. From our inversion procedure, we find that the V- to H-band scattering cross sections ratio increases 
^ ^ outside 40 AU, in line with the observed color gradient of the disk. This behaviour is expected if the grains have not been produced locally 
but placed in orbits of high eccentricity by a size-dependent pressure force, resulting in a paucity of large grains beyond the outer edge of the 
parent bodies disk. Because of the low luminosity of AU Mic, radiation pressure is inefficient to diffuse the smallest grains in the outer disk, 
even when the flares are taken into account. Conversely, we show that a standard, solar-like stellar wind generates a pressure force onto the 
dust particles that behaves much like a radiation pressure force. With an assumed M ^ 3 x 10^ Mq, the wind pressure overcomes the radiation 
pressure and this effect is enhanced by the stellar flares. This greatly contributes to populating the extended AUMic debris disk and explains 
the similarity between the fi Pictoris and AU Mic brightness profiles. In both cases, the color gradient beyond 120 AU for fi Pictoris and 35 AU 
for AUMic, is believed to be a direct consequence of the dust dynamics. 

Key words. Stars: circumstellar matter - Stars individual: AU Mic- Stars: flare ~ Planetary systems: formation - Scattering 

1. Introduction importantly the radial surface brightness profiles have similar 

shapes. Both dust disks are furthermore thought to be made of 

Recent visible and near-infrared coronagraphic imaging of the short-lived collisional debris and the gas to be largely dissi- 

nearby M-type star AU Mic revealed its geometrically and opti- p^jg^ (Liu et al. 2004b; Brandeker et al. 2004; Roberge et al. 

cally thin circumstellar dust disk (Kalas et al. 2004; Liu 2004a; 2005). AU Mic and ^6 Pictoris would moreover be co-eval since 

Krist et al. 2005; Metchev et al. 2005; Fitzgerald et al. 2005). jhey would belong to the same nearby stellar association aged 

The disk, seen almost perfecdy edge-on, extends over several 12+^ Myr (Barrado y Navascues et al. 1999; Zuckerman et al. 

hundreds of AU diameter The disk morphology presents inter- 20O 1 ) 
esting similarities with the well-studied disk around the young 

A5V main sequence star ^Pictoris: both disks are seen edge- Do these similarities between the two debris disks point to- 
on, imaging has evidenced radial and vertical asymmetries, and ward similar ongoing dynamical processes? The vertical asym- 

metrics of the /J Pictoris disk seen in scattered light have been 

Send offprint requests to: J.C. Augereau modeled assuming the presence of a planet on an inclined or- 
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Augereau et al. 2001). The AUMic disk appears sUghtly ver- 
tically distorted within 50 AU radius but current observations 
do not readily show large-scale vertical asymmetries simi- 
lar to those observed in the ySPictoris disk. Hypothetic plan- 
etary companion(s) around AUMic should therefore be on 
low inclination orbit(s) with respect to the disk midplane. The 
AU Mic disk also shows small-scale structures appearing as lo- 
cal brightness enhancements not present on each side of the 
disk (Liu 2004a; Krist et al. 2005). Such small-scale side-to- 
side asymmetries have also been reported forySPictoris (Pantin 
et al. 1997; Telesco et al. 2005) and could correspond to local 
enhancements of the collision rate due to planetesimals trapped 
with a planet (e.g. Roques et al. 1994; Telesco et al. 2005). 

More intriguing is the similarity between the midplane sur- 
face brightness profiles of the AU Mic and fi Pictoris disks in 
scattered light. The observed profiles both decrease almost lin- 
early with increasing the distance from the star then drop with 
power law indexes between -4.4 and -5.5. The break be- 
tween the two power-law regimes is observed around 35 AU 
for AUMic and 120 AU for /? Pictoris. According to Augereau 
et al. (2001), the surface brightness break around 120 AU in 
the y6 Pictoris disk corresponds to the outer edge of the colli- 
sional planetesimal ring that produces the observed dust par- 
ticles. The stellar radiation pressure of the central A5V star 
places the smallest bound dust particles on orbits of high ec- 
centricity with the same periastron as their parent bodies. These 
grains spend a large fraction of their time close to apoastron 
distributing this way the dust particles over a wide range of dis- 
tances far from the place they have been produced. This results 
in a disk that extends over hundreds of AU with a generic scat- 
tered light surface brightness profile further the 120 AU break 
and proportional to r"^ as observed (Lecavelier des Etangs et 
al. 1996; Augereau et al. 2001; Thebault & Augereau 2005). 
Because of the drag gas can induce on the dust particles, this 
model requires a low gas density to keep the scattered light 
profile consistent with the observations. Thebault & Augereau 
(2005) use this effect to place an upper limit on the gas density 
in the /? Pictoris disk. 

The previous model, successful to explain the scattered 
light brightness profile of the /? Pictoris disk, does not a pri- 
ori apply to the AUMic disk because of the low stellar lu- 
minosity of the central M-type star. The stellar radiation pres- 
sure on the grains is expected to be small and the steep slope 
of the scattered light profile further 35 AU would then reflect 
the profile of the parent bodies distribution rather than the ef- 
fect of the radiation pressure like for fi Pictoris (Kalas et al. 
2004). In this scenario, the similarities between the AU Mic and 
y6 Pictoris surface brightness profiles would just be coinciden- 
tal. This conclusion nevertheless depends on uncertain parame- 
ters. Radiation pressure depends on the grains properties (com- 
position, size) which are poorly constrained. AU Mic is also a 
flare star with UV and X-Ray rises that could periodically in- 
crease the radiation pressure force on the grains. Alternatively, 
a hypothetical stellar wind could create a corpuscular pressure 
force that would help difl'using the dust particles at large dis- 
tances. Finally, the collisional activity of the disk has not been 
precisely determined and requires the disk surface density to be 
properly estimated. 
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Fig. 1. Notations used in this study (see Sect. 2. 1) 

In this paper, we invert HST visible and ground-based near- 
infrared surface brightness observations to calculate surface 
density profiles. We derive in Sect. 2 the integral equation that 
links the observed profiles to the physical quantities we want 
to recover, and in Sect. 3 we detail our numerical approach. We 
then discuss the shape of the disk surface density as a function 
of the anisotropic scattering properties of the grains (Sect. 3.4). 
In Sect. 4 we constrain the grain size distribution based on the 
blue color of the disk in the visible (Krist et al. 2005) and we 
derive some disk properties (mass, optical thickness). Finally 
we argue in Sect. 5 and 6 that the observed r^^ surface bright- 
ness profile further 35 AU in scattered light is not coinciden- 
tal but instead reflects, like in the y6 Pictoris disk, the ongoing 
dynamics of the observed dust grains. Our scenario involves 
both the observed stellar flares and the presence of a wind drag 
force. 



2. Surface brightness profile 

2.1. Notations and assumptions 

We note (x, y, z) the observer's line of sight Cartesian frame 
with the star at the origin and with the x-axis pointing to the 
observer (Fig. 1). We consider a disk rotated around the y-axis 
by an angle / with / - corresponding to an edge-on disk. We 
assume -n/2 < i < 7r/2 to have cos/ > 0. We note (x',y',z') 
the Cartesian frame attached to the disk with the star at the 
origin and with the z'-axis perpendicular to the disk midplane. 
We let y' = y which implies that the y'-axis corresponds to the 
apparent major axis of the disk when / ^ 0. We furthermore de- 
fine (r, 6, z) and (r', 6',z') the cylindrical coordinates in the line 
of sight and disk frames, respectively. We let r^^^ the maximal 
radius of the disk and c/star the stellar heliocentric distance. 

The AU Mic disk is assumed to be optically thin in all di- 
rections and at all wavelengths which implies that the grains 
are directly illuminated by the central star This assumption is 
justified by the low fractional dust disk luminosity (Ldisk/i^* - 
6 X 10 '', Liu et al. 2004b), the low extinction of the star and the 
lack of central dark lane along the disk midplane in the visible 
and at near-lR wavelengths as regularly observed for optically 
thick edge-on disks. 
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To account for light anisotropic scattering by the grains, 
we introduce the scattering angle <p with <p - corresponding 
to forward scattering (see Fig. 1). Since r^^^ <K c/jtar^ we will 
assume in the following that - -cf) for the grains with z = 0. 



2.2. Scattered light surface brightness profile 

The surface brightness of an optically thin disk measured at a 
projected distance y from the star along the y-axis (i.e. at z = 0) 
writes 



B(y) 



CO 



p(/,6',z')^(r,(f))dx 



(1) 



where p(r', d',z') is the disk number density and <^(r, (f>) is the 
total flux measured at Earth emitted by a grain at position 
(r, 6--(f),z=0). This flux is a function of the observing wave- 
length but to simplify the notation, we will omit to indicate this 
dependence in the equations. 

At the distances probed by the visible and near-IR coro- 
nagraphic images, the thermal emission from the grains con- 
tributes little to the disk brightness and in the following we 
will only consider the single scattering of the stellar light by 
the dust particles. <t)(r, (p) reads then 0*crsca^"^ /(</>) where <£)* 
is stellar flux at Earth at the considered wavelength, osc-^ and 
f((/>) are respectively the scattering cross section and the scat- 
tering phase function of the grains averaged over the grain size 
distribution. Here we implicitly assumed that the line support- 
ing any incident stellar ray is an axis of symmetry for the light 
scattered by the grains (phase function depending only on (p). 
This is achieved for instance when the grains are spherical and 
homogeneous. 

We assume the disk is symmetric about the midplane and 
the y-axis. Equation (1) then re-writes 

. r°° , , , f(<p) + f(7i-(p) 

Jo '- 

= <5 o-sc^pir ,6 ,z) 5 TT'i'^ • (3) 

Jy r^ or' 

The infinite bound on r' in Eq. (3) can be replaced with 
the maximal disk radius r'^.^^ since the dust density is as- 
sumed to be null beyond r'^.^^. One can show that dxjdr' - 
r' cos"^ // "v/r^ - y^ for < < nil with 
y^/ cos /. We moreover have r^ = ( 
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• y^ sin I 
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cos^ / 



when z = 0. The scattered light surface brightness of the disk 
at a distance y from the star along the y-axis thus reads 



B{y) = O" 



t 



CT,,,p(r',0',z') 



f(c/>)+f(n-cf>: 

^ n -> ■ 2 ■ 



COSi 
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; r'dr' . 



(4) 



In the previous equation, we left the scattering cross section 
under the integral sign since the dust properties (composition 
and/or size distribution) could, in principle, depend on the dis- 
tance from the star. The grain size distribution in the /3 Pictoris 
disk is for instance expected to be distance-dependent (Fig. 
4 in Augereau et al. 2001) and the recent observations by 



Golimowski et al. (2006) support this model. The situation 
might be similar in the case of AU Mic as suggested by the 
observed gradual change of disk color in scattered light along 
the disk midplane (Krist et al. 2005; Metchev et al. 2005, see 
also Sect. 4.1). 



2.3. The case of edge-on geometrically thin disks 

The AU Mic disk is seen edge-on which simplifies Eq. (4) (/=0, 
r=r', 0=0'=-(p and z=z'=0). The brightness profiles shown in 
Fig. 2 have been obtained by integrating the disk surface bright- 
ness over a few pixels in direction orthogonal to the disk mid- 
plane. Therefore, B(y) as defined by Eq. (4) does not directly 
compare with the measured profile and the integration over z 
has to be taken into account. 

We can take advantage of the integration over z to get an es- 
timate of the surface number density of the disk instead of the 
number density in the disk midplane. The grains that contribute 
to S (y) in geometrically thin disks have their height z above 
the disk midplane that obeys z^ <s r^ which implies ^ -9. 
Because of the symmetries of the phase function we in addition 
have f(-6) = /(O) and f(7T + 6) = /(n-d) [see also Eq. (8)]. We 
approximate the scattered light surface brightness profile mea- 
sured along the disk midplane of an edge-on disk and integrated 
over Az/2 <K r in a direction perpendicular to the midplane by 



S(y)^0- 



r 

Jy 



(f(0) + f(n-0)) 



o-^cJ^jr, 0) 
r ^Jr^ - y^D,{r, 0) 



dr 



with S(r, 0) 
and Dj(r, 0) = Az 



p(r,0,z)dz 

oo 

S(r,0) 



/-Ma'^^'"'^'^)'!^ 



(5) 

(6) 
(7) 



where Y.{r,0) is the disk surface number density. Equation 5 
makes the link between the observed brightness profile S{y) 
and the disk properties (surface number density, phase func- 
tion and scattering cross section). It depends through D^i^r, 0) 
on the vertical height above the disk midplane over which the 
brightness profile has been integrated. Equation 5 belongs to 
the class of first-kind Volterra integral equations but because of 
the singularity m r - y, it actually belongs to the sub-class of 
Abel integral equations. 



3. Inversion of the AU Mic brightness profiles 

To infer the surface density profile of the AU Mic disk from ob- 
served scattered light profiles, we adopt a direct approach that 
consists in performing a basic numerical inversion of Eq. (5). 
This presents several advantages over a classical fitting ap- 
proach. In particular, this limits the number of free parame- 
ters since the inversion does not depend on any parametriza- 
tion of the surface density profiles. We can therefore effi- 
ciently explore the coupling between the surface density and 
the anisotropic scattering properties of the grains. 
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3.1. Observed surface brightness profiles 

In this study we consider two sets of scattered light bright- 
ness profiles iS{y)), both displayed in Fig. 2. The first set 
of profiles has been measured by Liu (2004a) in the H-band 
(Ac - 1.63 yum) at the Keck II 10-m telescope equipped with an 
adaptive optics system. A spatial resolution of 0.04", or 0.4 AU 
at the distance of AUMic (9.94 ± 0.13pc), has been achieved 
and the profiles have been derived from an aperture Az = 0.6" 
(~ 6 AU) wide in the direction perpendicular to the disk mid- 
plane. The Keck H-band coronagraphic images probe the disk 
down to ~ 15 AU. The second set of scattered light brightness 
profiles has been obtained in the visible with the HST/ACS in- 
strument by Krist et al. (2005). The ACS images have a reso- 
lution of ~ 0.63 AU and have been obtained with Az - 0.25" 
(~ 2.5 AU). The HST/ACS images present several features at 
distances closer than about 20 AU that are suspect according 
to Krist et al. (2005). The published ACS surface brightness 
profiles were multiplied by a factor of 1.125 to account for 
some confusion regarding the relation to the ACS geometric 
distortion and the photometric calibration points assumed by 
the HST Synphot package (Krist 2006, private communica- 
tion). 

The HST/ACS visible and Keck-AO near-IR images of the 
AUMic disk show some differences discussed in Krist et al. 
(2005). Their impact on the surface density profiles calculated 
with our inversion technique are addressed in Sect. 3.4. All ob- 
served brightness profiles were oversampled to reach a com- 
mon resolution (h in Appendix A) of 0.2 AU. The observed 
profiles were furthermore extrapolated with radial power laws 
beyond 65 AU for the H-band profiles and beyond 75 AU for 
the visible profiles (Fig. 2). The profiles in Fig. 2 have been ob- 
tained assuming O* = 12.7 Jy in the H-band and <C>* = 1.27 Jy 
in the HST/ACS V-band (F606W). 

3.2. Scattering phiase function 

The scattering phase function of solid particles can be very 
complex depending on their shape and composition. In order 
to efficiently explore the impact of the anisotropic scattering 
properties of the grains [term f(ff) + f{n - 6) in Eq. (5)] on 
the surface density profiles calculated with the inversion pro- 
cedure, we adopt the single parameter Henyey & Greenstein 
(1941) phase function 



/(0) 



l-g' 



4-n (1 + g^ - 2g cos cf)) 



3/2 



(8) 



where the parameter 8 - JL /(0) cos if> dQ is referred to as 
the mean asymmetry factor (g - for isotropic scattering and 
1^1 = 1 for purely forward/backward scattering). The ^parame- 
ter is related to the scattering properties of the individual grains 
through the following relation that ensures that the total scat- 
tered flux is preserved 



•J .Smi 



T^S g,scaX^dn(i)/crs, 



with cr^. 



•^■Smin 



rescadn(i) 



(9) 
(10) 













Distance in AU 

Fig. 2. Scattered light surface brightness profiles of the AUMic disk 
5(>') (divided by the stellar flux O* to aflow comparison) in the H- 
band (Liu 2004a) and in the ACS V-band (Krist et al. 2005). Beyond 
r = 65 AU, the original H-band profiles have been extrapolated with 
radial power laws (NW: r^*-^, SE: r"* '). The visible profiles have been 
similarly extrapolated beyond 75 AU (NW: r""*^, SE: r"''). The ver- 
tical shift between the two profiles has two origins: the disk color in 
scattered light (Sect. 4.1), and the difi"erence in spatial extent (Az) in 
a direction orthogonal to the disk midplane used to measure the disk 
brightness. 



where s is the grain size, s^nin and i^ax are the minimum and 
maximum grain sizes, respectively, and dn(s) is the unity nor- 
malized differential grain size distribution (J dn(s) - 1). 
gsca is the dimensionless scattering efficiency and ^ is the 
asymmetry factor that depends on the grain properties (size, 
composition, structure, ...) and on the wavelength. 

To efficiently explore the impact of scattering anisotropy on 
the results and to limit the number of free parameters, we will 
assume that the g parameter is a constant. The inversion algo- 
rithm can nevertheless account for mean asymmetry parame- 
ters that depend on the grain position in the disk. We will also 
only consider the absolute value of g since we cannot distin- 
guish between positive and negative values due to the edge-on 
orientation. The asymmetry factors \g^\ and |g^| in the H-band 
and ACS-V-band, respectively, are hence the only free physical 
parameters when reconstructing surface density profiles from 
scattered light profiles using Eq. (5). 

3.3. Vertical and azimuthal profiles 

Following Krist et al. (2005), we assume that the vertical profile 
of the disk density can be approximated by a Lorentzian func- 
tion (l + z^/H^(r)] . The scale height H(r) is parametrized as 
foflows: //(20AU) = 0.8675 AU, H(r) oc /""* for r < 49 AU, 
H(r) oc r^^ for 49 < r < 80 AU and H(r) oc r"^ for r > 80 AU. 
D^(r, 0) thus writes 



D,(r,fl) = Az^arctan-'(^ 



such that D,(r, 0) = Az for Az » 2 H(r). 



(11) 
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Fig. 3. Surface density profiles obtained by inverting the Liu (2004a) H-band observations and the Krist et al. (2005) HST/ACS visible observa- 
tions of the AU Mic disk (top and bottom panels, respectively, for each figure), assuming four mean scattering asymmetry parameters. The solid 
line corresponds to the NW side; the dotted line to the SE side. These results have been obtained assuming ^ = 1 (see Sect. 3.3 and App. B) but 
the profiles are poorly sensitive to the assumed azimuthal profile provided 6 < 5. 



The scattered light images of the AU Mic disk evidence 
side-to-side asymmetries. Also, to appreciate how the inver- 
sion of the NW brightness profile affects the inversion of the 
SE brightness profile (and reciprocally), we introduce an ad 
hoc dependence on 6 in the inversion process and we pro- 
ceed to the inversion of the NW and SE brightness profiles si- 
multaneously. We re-write the dust surface number density as 
S(r, ff) = So(r)0(r, 6) where 0(r, ff) is a unit-less function given 
in Appendix B and characterized by a single, purely geometri- 
cal, parameter 6. The role of the 0(r, ff) function is to smoothly 
link the NW and SE sides of the disk. During the inversion 
process, the surface density is calculated approximately along 



the y-axis (see Fig. 1). The 0(r, ff) function forces the NW and 
SE surface density profiles to be identical along the x-axis and 
equal to their mean value. 

3.4. Surface density profiles 

To reconstruct the surface density profiles of the disk, we have 
used a classical numerical scheme to invert Eq. (5). The algo- 
rithm is detailed in Appendix A in the case of edge-on disks 
but it can easily be extended to the case of inclined disks (e.g. 
the ring about the F5 star HD 181327, Schneider et al. 2006). 
As a sanity check, the calculated surface density profiles were 
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systematically incorporated into our 3D debris disk model to 
compute synthetic brightness profiles (Augereau et al. 1999, 
2001). These synthetic profiles were always found to be in ex- 
cellent agreement with the original observed profiles. 

The results of our inversion procedure are reported in 
Figure 3. This figure displays the various surface density pro- 
files found assuming four different values for the mean asym- 
metry parameter \g\. For each \g\ value, the top panel shows 
the surface density profile that reproduces the H-band scattered 
light surface brightness profile measured by Liu (2004a), while 
the profile in the bottom panel is derived from the ACS-V-band 
observations (Krist et al. 2005). The surface density profiles are 
found to depend quite sensibly on the assumed \g\ value. Large 
\g\ values promote the emission at scattering angles close to 
or TT (depending on the sign of g) which is achieved when the 
grains are at distances from the star large compared to their 
projected distance (r » y in Eq. (5)). The amount of dust at 
large distances thus tends to increase with increasing |^| values 
and the regions interior to the disk thus become more and more 
dust-depleted. As a result, the disk as derived from the H-band 
brightness profiles resembles a ring peaked around 35 AU for 
sufficiently large \g^\ values. Furthermore, |^| cannot be signifi- 
cantly larger than ~ 0.6 since for larger |^| values, the inversion 
algorithm requires the surface density profiles to be negative at 
some radii to compensate for the too large flux, compared to 
the observations, intercepted on the line of sight. 

There are obvious diflFerences between the surface density 
profiles obtained from the two data sets used in this study. They 
basically reflect the diflFerences between the original images 
(see discussion in Krist et al. 2005). Along a given arm, the 
same density features can generally be identified (though often 
shifted by 1-2 AU with respect to each other) but not with the 
same strength. For instance, the NW surface density profiles 
inferred from the H-band observations, peak at ~ 37 AU and 
show a smaller peak around 30 AU, while the feature strengths 
are reverted in the case of the surface density profiles obtained 
from the visible images. These two peaks correspond to fea- 
tures (D) and (C), respectively, according to the nomenclature 
of Liu (2004a), and feature (D) is only marginally resolved 
in the ACS images. Along the SE side, features (A) and (B) 
around 25-27 AU and 33-34 AU, respectively, show strength 
ratios that depend on |^| and on the data set used to calculate 
the disk surface density. The visible images in addition har- 
bor an extended clump around 50 AU that results in a broad 
feature at that distance, and which is only marginally present 
in the surface density profiles deduced from the H-band ob- 
servations. Almost all profiles, nevertheless, show a common 
dust depletion around 22-24 AU, the exact position depending 
on the data set used to calculate the surface density profiles. 
Besides the differences, we note that all profiles show complex 
structures that can hardly be approximated with simple power 
laws. 

The differences between the surface density profiles de- 
rived from the visible and infrared observations may result 
from a combination of instrumental effects and grains optical 
properties varying throughout the disk midplane. The inver- 
sion of additional observed profiles and more constraints on 
the grain properties would help to disentangle between these 
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Fig. 4. Scattering cross section ratio cr^a/''"sca obtained by dividing 
ttie profiles calculated for |g^| = |g^| = 0.4 and shown in Fig. 3. The 
gray area corresponds to the range of scattering cross section ratios 
inferred from the analysis of the disk color in the ACS bands (Sect. 4. 1 
and Table 1). The R-H colors inferred by Metchev et al. (2005) are 
displayed as two black dots at 20 and 55 AU. 



two possibilities. But if we assume that the differences are real, 
then the ratio cr^^Jlo(r) / crl^^Jloir) is an estimate of the V- to H- 
band scattering cross section ratio as a function of the distance. 
This ratio is plotted in Figure 4 in the case \g^\ - |g^| - 0.4. 
This suggests that the cr|^a/cr^a ratio gradually increases be- 
yond ~ 40 AU, in line with the color gradient of the disk over 
the same region (Krist et al. 2005). Qualitatively, this would 
indicate that the large grains become less and less abundant, 
relative to the smallest grains at the same position, with in- 
creasing the distance from the star. This behaviour is expected 
if the grains have not been produced locally but placed in orbits 
of high eccentricity by a pressure force oc s^ jiP' + z^) (where s 
is the grain size) as shown in Fig. 4 of Augereau et al. (2001). 
We discuss in Sect. 5 and 6 the origin of pressure forces in the 
case of AU Mic. 



4. Grain and Disk Properties 

4.1. Disk color and grain size distribution 

The output of the inversion procedure is not the actual surface 
density profile of the disk, 'Loir), but its product by the mean 
scattering cross section osc-^. To evaluate cr^^;^, and hence Eo(r), 
we search for the grain properties that can reproduce the blue 
color of the AU Mic disk reported by Krist et al. (2005) in the 
HST/ACS bands. These authors find F435W/F606W flux ra- 
tios, relative to the star and averaged over a 43 AU long aper- 
ture centered at r = 42AU, of 1.13 ± 3% and 1.15 + 3% for 
the NW and the SE sides of the disk, respectively. The mea- 
sured F435W/F814W flux ratios amounts to 1.35 + 3% (NW) 
and 1.44 ± 3% (SE). Qualitatively, the blue color indicates that 
the visible scattered light images are dominated by elementary 
(sub-)grains of the order or smaller then a few tenths of mi- 
crometers. 
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Fig. 5. Minimum grain sizes .s„,in and power-law indexes k of the size distribution consistent with the HST/ACS colors measured by Krist et al. 
(2005). Gray area: F435W/F606W, dark area: F435W/F814W, dashed area: intersection. The cross sign corresponds to the nominal pair (iViim 
k) that produces a perfect match to the two colors. The nominal values are given in Table 1. 

Table 1. Grain properties inferred from the interpretation of the disk colors in the visible (Krist et al. 2005) and implications for the disk mass. 





Parameter*"* 


Astronomical Silicates 


Graphite Grains 




Description 


NW 


SE 


NW 


SE 


[unit] 


Nominal minimum grain size 
Nominal power law index 
H-band mean asymmetry factor 
V-band mean asymmetry factor 
H-band mean scattering cross section 
Ratio of mean cross sections 
V-band mean albedo 


K 

lg"l 

\s'\ 


0.115 

-3.74 
0.637 
0.675 
16.0 
2.05 
0.795 


0.123 
-4.01 
0.604 
0.661 

12.7 
2.56 
0.822 


0.0583 
-3.95 
0.358 
0.455 
2.12 
2.15 
0.518 


0.0609 
-4.18 
0.318 
0.421 
1.64 
2.65 
0.513 


[jum] 
[10-'" cm^] 




5 "a 

Mdisk 




H-band data set 






Disk scattering cross section*' 
Dust mass*"^' 


4.93 - 6.97 
2.60-3.67 


4.51-6.30 
0.952- 1.33 


4.93 - 6.97 
0.583 - 0.824 


4.51-6.30 
0.306 - 0.428 


[AU2] 
[10-^Me] 


Total dust disk Mass'*^* 


3.55- 


-5.00 


0.889 - 


- 1.25 


[lO-^'Me] 




Mdisk 




V-band data set 






Disk scattering cross section*' 
Dust mass*'^' 


15.9-23.7 
4.10-6.10 


15.5-23.1 
1.28-1.91 


15.9-23.7 
0.876-1.31 


15.5-23.1 
0.398 - 0.593 


[AU^] 
[10-'*Me] 


Total dust disk Mass**^* 


5.38- 


-8.01 


1.27- 


1.90 


[10-^ Me] 



" H and V superscripts stand for H-band and ACS V-band, respectively 

'' for < Igl < 0.6 and S = I. Sg^a, given by Eq. (15), does not depend on the assumed grain composition 

° mass inferred from the analysis of the scattered light observations. See Sect. 4.3 for an estimate of the dust mass deduced from the analysis 

of the disk thermal emission. Mass obtained assuming s^^ax = 1 nrni, psiUcaies = 3.5 g.cm"^ and Pgrapwte = 2.2 g.cm"^. The total dust disk 

mass is the sum of the NW and SE dust masses. 



Provided the phase functions are not too different, the mea- 
sured colors are only dependent upon the ratios of mean scat- 
tering cross sections in the ACS passbands. We have consid- 
ered optical constants for astronomical silicates (Weingartner 
& Draine 2001) and graphite grains (Laor & Draine 1993), and 
we used the Mie theory, appropriate for hard spheres, to calcu- 
late the scattering efficiency Qsca- We have assumed a differen- 
tial grain size distribution of the form dti(s) oc s''ds with k < 
between imin and imax. The maximum grain size is assumed to 



be large enough to not affect the results. The scattering cross 
sections were calculated using Eq. (10). 

We have searched for the pairs (imin, k) that give scattering 
cross section ratios equal to the two HST/ACS colors. Each side 
of the disk has been considered independently. The results are 
reported in Table 1 and in Fig. 5, where the cross sign indicates 
the nominal couple (Smin^ k), namely the couple that produces 
a perfect match to the two ACS disk colors. For astronomical 
silicates, the blue color of the disk is reproduced when the mini- 
mum grain size Sniin is smaller than 0. 14jum with a nominal imin 



J.-C. Augereau & H. Beust: The AU Mic debris disk 



value of 0.12/im. For graphite grains, we find «„;„ < 0.07 /im 
with a nominal s^i„ ^ 0.06 ;um. For a given grain composi- 
tion, our simple model predicts quite similar nominal s^^m val- 
ues for both the SE and NW sides of the disk, but the nominal 
slopes K differ: the size distribution in the NW side is found to 
be less steep than SE size distribution where the disk is found 
to be even bluer (Krist et al. 2005). k values consistent with the 
HST/ACS colors are generally found to be slightly smaller than 
the canonical -3.5 value that holds for systems in collisional 
equilibrium (Dohnanyi 1969). The minimum grain sizes we in- 
fer for astronomical silicates are smaller than the s^in value 
of 0.5;^Q2yum obtained by Metchev et al. (2005) based on the 
R-H color of the disk measured around 50-60 AU and assum- 
ing a- = -3.5. 

The mean asymmetry factors reported in Table 1 were cal- 
culated using Eq. (9). For astronomical silicates, \g\ reaches or 
slightly overcomes the 0.6 limit derived from the inversion 
of the scattered light profiles (Sect. 3.4). For more refractory 
grains such as graphite grains, the mean asymmetry factor re- 
mains smaller than the \g\ = 0.6 limit and is close, in the V- 
band, to the 0.4 value inferred by Krist et al. (2005) from the 
disk modeling. 
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Fig. 6. Vertical optical thickness assuming astronomical silicates (see 
Table l).Top panels: surface density profiles inferred from the inver- 
sion of the H-band observations; bottom panels: surface density pro- 
files derived from V-band observations. 



4.2. Disk optical tliickness 

In the previous sections, we have calculated the surface density 
and grain properties of the AU Mic disk assuming an optically 
thin medium at all wavelengths and in all directions. This as- 
sumption has been justified in Sect. 2. 1 by the low fractional 
disk luminosity and the appearance of the disk at visible and 
near-IR wavelengths. We verify in the following that our re- 
sults are self-consistent. 

The vertical optical thickness of the AU Mic disk, Tj_(r), 
reads 

^o(r) = ^-^ X ^^,,So(r) (12) 



1 



where co^ is the albedo of the dust particles in the V-band av- 
eraged over the grain size distribution. The ratio cr|^^/cr^^ and 
the mean albedo w^ have been estimated assuming the grain 
properties summarized in Table 1 . In the visible, the disk ver- 
tical thickness T^(r) does not exceed ~ 5 x 10"^ indicating 
that the disk is vertically optically thin (Fig. 6). We checked 
that with the assumed grain properties, the vertical thickness is 
maximum in the visible implying Tj_ <K 1 at all wavelengths. 

The midplane density of the disk writes Eo(r)/(7r//(r)) and 
the optical thickness along the disk midplane thus reads 
^v 



T||('") 






(13) 



TH(r) J,„, nH(r) 

where rmin = 15 AU in the H-band and r^in = 20 AU in the 
ACS V-band (see Figs. 2 and 3). Figure 7 shows that T||(r) does 
not exceed ~ 4 x 10"^ in the visible whatever \g\ < 0.6 and 
this result is independent of the data set used to calculate the 
surface density profiles. Given that t± and thus Ty both reach 
their maxima at visible wavelengths in this model, these results 
indicate that the disk is optically thin in all directions and at all 
wavelengths. 
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Fig. 7. Midplane optical thickness assuming the same grain properties 
as in Fig. 6. 



4.3. Disk mass and implication for the grain properties 

From the surface density profiles calculated in Sect. 3.4 and the 
constraints on the grain size distribution obtained in Sect. 4.1, 
we can estimate the mass of the dust disk Mdisk further out than 

^min 



_ ^sca(lgl) ^, 
Mdisk - X 






5 dn(5) 



(14) 



where p is the grain mass density. The disk scattering cross 
section 5sca(l^l) writes 



5sca(l§l) 



I ' 



crsca2(r, G)rArAd 



(15) 



J.-C. Augereau & H. Beust: The AU Mic debris disk 



where r^^x was fixed to 140 AU. The results are given in 
Table 1. 

The total mass of dust derived from the scattered light anal- 
ysis is of the order, to a few times, the mass of the Ceres as- 
teroid (~ 1.46 X 10"'* Me) in the Solar System. This is about 
10 to 80 times the estimated dust mass in the Kuiper Belt (~ 
10"^ Me, Backman et al. 1995). But it is significantly smaller 
than the dust mass calculated by Liu et al. (2004b) from sub- 
millimeter observations (1.1 X 10"^ Me). In fact, a dust mass as 
low as a few lO"'* Me is not expected at an age of about 10 Myr 
(e.g. Najita & Williams 2005). 

To help identifying the origin of the mass diff'erence, we 
calculated the disk Spectral Energy Distribution (SED) assum- 
ing the surface density profiles and the grain properties ob- 
tained from the analysis of the scattered light observations. We 
assumed the star is in a quiescent period (no UV-flares, see 
Fig. 9 and Sect. 5.2). We find that the total 25 /im flux is al- 
ways well reproduced while the 60-70 //m fluxes are system- 
atically underestimated by factors of 1 .5 to 2. More critically, 
the observed sub-millimeter emissions are 10 to 20 times larger 
than predicted when silicates are considered. These flux ratios 
amount to 15 to 30 for the graphite grains. 

On the other hand, we find that the whole disk SED can be 
reproduced provided the constraints on the grain size distribu- 
tion are relaxed. Good fits to the SED are for instance obtained 
for Sjnin - 1.0 /im and k ^ -3.4 when silicate grains and the 
surface density profiles derived from the H-band images are 
assumed (Fig. 9). With the V-band surface density profiles, we 
find Smin - 0.1 //m and k ^ -3.4. In both cases, the dust disk 
mass is approximately 7 x lO""' Me, in better agreement with 
the submillimeter dust mass (1.1 x 1 0"^ Me) calculated by Liu 
et al. (2004b) assuming a single temperature formula. 

The grain size distributions that reproduce the disk colors in 
scattered light do then underestimate the amount of large grains 
in the disk, which explains the disagreement between the two 
independent mass estimates. It might be that the dust properties 
in Table 1 are only valid in a small range of grain sizes. This 
could be the case, for instance, if the actual grain size distribu- 
tion departs from a simple power law as expected for systems in 
collisional equilibrium (see discussion in Thebault et al. 2003). 
It might also be possible that the grains have complex struc- 
tures, and hence optical properties, that cannot be consistently 
reproduced over the whole spectrum assuming hard spheres. 



5. Diffusion of dust particles in tlie outer disk. 
I. Time-variable radiation forces 

Radiation pressure is a key parameter controlling the dynamics 
of dust particles in usual debris disks. The grains are thought to 
be produced by numerous, colliding or evaporating planetesi- 
mals orbiting the star on nearly circular orbits. Depending on 
their size, the radiation pressure can transport the grains far 
away from their production zone and fill the outer regions of 
the disk. This model explains the observed radial profiles of 
the yS Pictoris and HD 181327 disk for instance (Lecavelier des 
Etangs et al. 1996; Augereau et al. 2001; Thebault & Augereau 
2005; Schneider et al. 2006). We discuss below and in the next 



Section whether this model could or not apply to the AU Mic 
disk. 



5.1. Equation of motion and dynamical time-scales 

A dust particle of mass m and geometric cross section cr feels 
a full force from the stellar radiation field of energy flux F that 
reads 



dV 



Fo-Qp, 



(■-9-; 



(16) 



(Burns et al. 1979), where r is the position vector of the par- 
ticle with respect to the star, v its velocity vector, m = r/r is 
the radial unit vector, r = ii ■ v is the radial velocity, and c 
is the speed of light; Qp^ is the dimensionless radiation pres- 
sure efficiency coefficient averaged over the stellar spectrum. 
In this expression, the velocity independent term (Fo-Qpr/c) u 
is often referred to as the "radiation pressure", while the other 
terms constitute the Poynting-Robertson drag (hereafter PRD) 
(Gustafson 1994). The radiation pressure is radial, and it is usu- 
ally described by its ratio Pp^ to stellar gravity F^ - GMm/r^, 
where M is the mass of the star and G the gravitational con- 
stant. As F cc r"^, /3pr does not depend on r for a given particle. 
Equation (16) is then rewritten as 



dV 
d? 



jSprGM 



PR 



('-;)-; 



(17) 



As V <s c, the radiation pressure largely overcomes the PRD. 

The long-term effect of PRD is to cause a decrease of the 
orbit. Computing the rate of change of the specific angular mo- 
mentum due to PRD, and assuming that the orbit is circular, we 
derive the orbital decrease as 

dr iPpfiM a 

-T-^-— ^-2-, (18) 

df re r 

where a is a constant parameter. Numerically, we have 

Q'=/3prX6.24xlO"'*— -AU^yr"', (19) 

a value already quoted by Wyatt (2005), and the characteristic 
PRD time-scale writes 

fpRD = k/(dr/dOI=0.5rV- 

As we see below, /3pr barely reaches 0.2 for submicron-sized 
particles. Assuming ySpr = 0.1 and M = 0.59 M© (Houdebine 
& Doyle 1994), the PRD time-scale turns out to range between 
2.7 X 10'' yr at r = 20 AU and 7 x 10^ yr at 100 AU. As quoted 
by Kalas et al. (2004), at sufficiently large distance it is longer 
than the age of the system (12^^ Myr, Barrado y Navascues et 
al. 1999; Zuckerman et al. 2001). 

The PRD time-scale moreover needs to be compared to the 
collisional time-scale, as collisions between grains of various 
sizes will tend to prevent the orbital decrease of those for which 
PRD is more efficient. The collision time-scale may be esti- 
mated as 



P/2 



fcoU 



^{s^) 2(r,,j:o(r)^/GM 



(20) 
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Fig. 8. Collisional time-scale assuming the surface densities displayed 
in Fig. 3 and the grain size distribution for astronomical silicates from 
Table 1. 



where 7r(i^) is the mean geometric cross-section of the grains 
averaged over the grain size distribution (Backman & Paresce 
1993). The result is shown on Fig. 8 for astronomical sili- 
cates. We see that for any distance r, the collision time-scale is 
more than one order of magnitude less that the PRD time-scale. 
Hence we may safely stress that collisions dominate over PRD; 
if any effect of PRD is to be investigated, it must be averaged 
over the size distribution, i.e. assuming an average ySpr ^ 0.01. 
The corresponding PRD time-scale turns out to be larger than 
the age of the star, so that we may assume that PRD is of negli- 
gible importance in the dynamics of the AU Mic disk. Actually, 
the same conclusion is derived by Wyatt (2005) for all the de- 
bris disks known today. 

Assuming the Poynting-Robertson drag force on the grains 
is negligible in the AU Mic disk, the equation of motion of a 
grain thus reads 



dV_ GM(l-;6pr) 
df2 r^ 



(21) 



This is a classical equation of a Keplerian motion. While suf- 
fering a radiation pressure force, a particle keeps following a 
Keplerian orbit around the star, but with an effective central 
mass M(l -jSpr). Depending on the value of jSpr, that orbit may 
be very different from the original orbit of the parent body. It is 
eccentric, shares the same periastron as the parent body orbit, 
but has a very different apoastron. At release time, the grain as- 
sumes the same orbital velocity as its parent body. If the orbit 
of the parent body is circular, energy and angular momentum 
preservation impose the grain orbit to obey 



a - 



a(l -ygpr) 
1 - 2;6pr 



/^P 



l-/?pi 



e' = 



l-2/3p 



(22) 



Here a is the semi-major axis of the original orbit, and a' is that 
of the resulting dust particle orbit; e' is the eccentricity of that 
orbit and Q is its apoastron. If ySpr is sufficiently close to 0.5, 



the grain is transported at a large distance thanks to a large Q 
value. If j6pr > 0.5 it is expelled from the system. 



5.2. The spectrum of a AUMic, a flare star 



The ySpr ratio depends on the grain properties (size, composi- 
tion) and on the SED of the central star. For any given star, /3pr 
is maximum for grains with radius s in the range 0.1-1 fim, falls 
off like i"' (geometric regime, gpi- = constant) at larger sizes, 
and tends to reach a constant for smaller grain sizes (Q^^ ^ s, 
Rayleigh scattering regime, Gustafson 1994). In the solar case, 
ySpr barely reaches 1 for submicron-sized particles (Gustafson 
1994; Beust et al. 2001), but in the case of classical Vega-like 
stars (A-type stars), j3pr reaches high enough values to enable 
ejection (a few micrometers in the case of jSPictoris depend- 
ing somehow on the grain composition, e.g. Artymowicz 1988; 
Beust et al. 2001; Augereau et al. 2001). 

Contrary to those stars, AUMic is an M-type star; hence 
radiation pressure is expected to be weak (Kalas et al. 2004; 
Liu 2004a; Krist et al. 2005). In order to compute it, we need 
to model the SED of the star. We fit the stellar spectrum 
with a NextGen model (Hauschildt et al. 1999) with \og(g) - 
4.5 (CGS) and an effective temperature of 3700 K. This gives 
a stellar luminosity of 0.092 L©. However, this model is not 
enough for describing the full SED of AU Mic. 

AUMic is known to frequently present X-ray and EUV 
flares (Kunkel 1973; Cully et al. 1993; Tsikoudi & Kellett 
2000; Magee et al. 2003). The flaring rate is estimated by 
Kunkel (1973) to 0.9 per hour. During these flares, the EUV 
and X-ray luminosities increase by a typical factor ~ 10 (Cully 
et al. 1993; Monsignori-Fossi et al. 1996; Tsikoudi & Kellett 
2000). We may estimate that the star spends ~ 10 % of the time 
in the flaring state. Even in quiescent state, the coronal activ- 
ity of the star is not negligible. Tsikoudi & Kellett (2000) note 
for the X-ray and EUV luminosities Lx = 2.24 x 10^^ ergs"' 
and Leuv = 2.9 x 10^*^ erg s"^ Hunsch et al. (1999) gives 
Lx = 5.5 X 10^^ ergs"'. Mitra-Kraev et al. (2005) reports si- 
multaneous X-ray and UV observations of AU Mic with XMM- 
Newton and give luminosities in both domains. Their X-ray lu- 
minosity is compatible with the previous estimate of Tsikoudi 
& Kellett (2000). These values are far above what should be 
expected from the stellar model described above. The X-ray 
and EUV luminosities may contribute to the radiation pres- 
sure felt by the dust particles, so that it is necessary to take 
them into account in the stellar SED model we use. The global 
shape of the X-ray and EUV spectrum of AU Mic is given in 
Monsignori-Fossi & Landini (1994); in the UV spectral domain 
between 1 150 A and 2000 A, lUE data are available (Landsman 
& Simon 1993); Archive HST/FOS spectra cover the spectral 
domain between 1600 A and 6800 A. Finally, broadband data 
(U, B, V, Re, Ic, 2MASS, IRAS, Spitzer, SHARC II, SCUBA) 
provide additional measurements towards the infrared (Chen et 
al. 2005, and references therein). 
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Fig. 9. Observed SED of the AU Mic system. The solid black curve corresponds to the NextGen stellar atmosphere spectrum that fits the visible 
and near-infrared data. The dashed and dashed-dotted curves are the SED fits for the quiescent and flare states respectively (Eq. 23). Two fit 
examples to the disk SED are also overplotted assuming the two surface density profiles obtained for |g| = 0.4 (see Sect. 4.3) and silicate grains. 



Given these constraints, we fit the full spectrum of AU Mic 
by adding to the NextGen stellar atmosphere model the follow- 
ing flux, in Jansky at Earth 



F(A) = exp 



with 



fl() + 



as 



1 + 



ln(/l) - a I 

02 



^3 (22^4 



2i2 



(^ in A) ,(23)==- 



flO = 


-- -13.7728, 


fli 


= 8.9395, fl2 = 1.3165, 


03 -- 


-- 39. 9113 K 


fl4 


= 12.5702 A, fl5 = 13.2092 



(24) 



This function fits accurately the various data in the UV and 
EUV domains and is compatible with the X-ray and EUV lumi- 
nosities given by Tsikoudi & Kellett (2000). This holds for the 
quiescent state. During flares, we expect the flux in X-ray and 
UV to be higher. There is nevertheless good evidence for syn- 
chronicity between flare events in X-ray and in UV (Tsikoudi 
& KeUett 2000; Mitra-Kraev et al. 2005), the time-lag between 
the two spectral domains being less than a few minutes (Mitra- 
Kraev et al. 2005). Hence we decided to mimic flare spectra by 
multiplying the F(A) function we add to the NextGen model 
by a fixed, average factor we estimate to 10 (see Tsikoudi & 
KeUett 2000), for wavelengths shorter than 2500 A'. The re- 
sult is displayed in Fig. 9, where we superimpose the fits to the 
NextGen model and to the various measurements. We see that 
shortwards to ~ 0.3 jum, the resulting spectrum is significantly 
above the NextGen model and will have an influence on thejSpr 
parameter felt by the various grains. 



' in practice, we multiplied Eq. (23) by 10 x (0.5 - arctan((/l - 
2500)/500)/;r) with A in A 
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Fig. 10. /3p, ratio for astronomical silicates (black curves) and graphite 
grains (grey curves) assumed to be spherical and homogeneous, as a 
function of their radius. The dashed curves correspond to the quiescent 
state of the star, the dotted curves to the flare state, and the sofid curves 
to the temporal average. 



5.3. Radiation pressure during X-ray and UV flares 

As the X-ray and UV spectrum of AUMic is subject to fluc- 
tuations due to flares, the /3pi- ratio for a given grain may be 
time- variable. Hence the dynamics of the dust particles may 
be more complex than described above where we assumed a 
constant /3pi- value. As the flares are very frequent, we cannot 
neglect this effect. We thus investigate now the dynamics of 
solid particles with a time- variable ySpi. 

In a simple description, we assume that the star periodically 
alternates quiescent and flare phases, respectively characterized 
for a given grain by ^Spi-q and ySprj ratios, with ySpr,f > ySpr,q. 
The star is assumed to spend a fixed fraction y of the time 
in flare phases. Strictly speaking, the dust particle will follow 
bits of Keplerian orbits alternatively characterized by dynami- 



12 



J.-C. Augereau & H. Beust: The AU Mic debris disk 



cal masses M(l -j6pr,q) and M(l -ySprf). It can be shown (see 
appendix C) that the global effect is the same as if fi^^ was 
constant, equal to its temporal mean. The dust particle keeps 
following a Keplerian orbit around the star, but characterized 
by a y6pr,eff ratio equal to the temporal average of ySpr, namely 

y6pr,eif = (1 - r)/?pi-,q + T/^pr.f- 

Our problem reduces therefore to estimating the ySpr.eff ra- 
tio, i.e. the parameters y6pr,q, y6pr,f and y. Figure 10 shows the 
value of ySpr as a function of the grain radius, computed from 
the SED of Fig. 9, in both quiescent and flare phases, for two 
grain compositions and assuming hard spheres. We also plot 
thCySpr^efF average, assuming y - 0.1. Figure 10 shows that un- 
less the grains are very refractory, they do not sufl'er in qui- 
escent phases a radiation pressure force large enough like in 
the p Pictoris disk to a priori produce an extended dust disk. 
For silicate grains, ySp, barely reaches 0.1 for submicron-sized 
grains. The radiation pressure is then only able to transport the 
grains up to 25% of their original distance at most (Q' - 1.25a 
forySpr = 0.1). In order to reach a distance of 210 AU (maximal 
disk radius, Kalas et al. 2004), a grain produced at 40 AU needs 
ySpr - 0.405, a value never reached during quiescent phases 
even for refractories. During flares conversely, the ySpr ratio of 
graphite grains can reach 0.5; but the temporal meanySpreff does 
not. With 7 = 0.1 (estimated from the light curves of Tsikoudi 
& Keflett 2000; Magee et al. 2003; Mitra-Kraev et al. 2005), 
PptfiS remains close to y6pr,q, and the conclusion derived above 
during quiescent phases still holds. The flares help diffusing the 
particles slightly further out, but not enough for explaining the 
brightness profile of the disk beyond ~ 35 AU. Another expla- 
nation must be found. 



6. Diffusion of dust particles in the outer disk. 
II. Wind Pressure Force 

Another possibility for blowing dust particles into the outer 
disk is the interaction with a stellar wind. Since AU Mic is a 
late type star, a solar-like coronal wind is likely to arise from 
the star. This is supported by the coronal activity of AU Mic. 



6.1. Dynamics of grains in tine presence of a wind 

The force due to a flux of protons arising from the star is 
very similar to the radiation force, and may be described as 
(Gustafson 1994) 



(1994). We write 



dV 



sw 






(''^y-l 



(25) 



where Op is the momentum density of the proton flux at the 
considered stellar distance, and Co is the free molecular drag 
coefficient. This equation is of course very similar to Eq. (16), 
and as for the radiation pressure case, the velocity indepen- 
dent term (i>pO-(Cu/2)u is dominant. Logically, we define the 
dimensionless ratio ySjw of that force to stellar gravity, so that 
the equation of motion becomes now similar to Eq. (17) with 
jSpr changed to ySsw The ySsw ratio has been introduced in the 
solar wind case by Mukai & Yamamoto (1982) and Gustafson 



Aw = 



%(t(Cd/2) 
GMm/r^ 



(26) 



The main difference between ySpr andySsw is that the latter is not 
necessarily independent from the stellar distance r. We inves- 
tigate this dependence. If psw and Vsw are the mass density and 
velocity of the wind at distance r, we have Op = psw^s,,. The 
coefficient Co is given by Gustafson (1994): 

1 



Cd = 



25^+1 ,2 45^-1-45^ 



2(1 -e)V^ 
"^ 3S~ 



2S^ 



erf (5) 



&(i + y.^ 



(27) 



where S = ^jm^l{2kT^)u is the ratio of the relative veloc- 
ity M to the thermal velocity of the protons -sJlhTp/rrip (k 
is Boltzmann's constant), e is the fraction of re-emitted pro- 
tons at dust temperature T^, and Fp is a coefficient intro- 
duced by Mukai & Yamamoto (1982) to account for sputter- 
ing of molecules of mass nis. In practice, for any standard wind 
model, at sufficiently large distance we have u ^ Vs„ (we ne- 
glect the orbital velocity), and for S » 1, Co is close to 2. 

If we assume spherical symmetry, then we may write the 
continuity equation M = 4;7rr^pswV'sw where M is the mass loss 
rate. If we assume that the grain is spherical, with radius s and 
density p, then the ratio ySsw may be rewritten as 



Aw = 



3 Mv.wCp 
327r GMps 



(28) 



Thus ySsw may be a function of the distance to the star as Vws 
and Cd are not necessarily constant throughout the disk. The 
stellar wind blow-out size then reads (,Ss„ = 0.5) 

3 MvswCd . . _ - ,-„, 

isw = 7^-^777 — withCD-2. (29) 

167r GMp 



6.2. Wind model and Mass loss rate 

If we want to describe further we need to adopt a wind model. 
The most standard model for a solar-like coronal wind is the 
transonic isothermal solution by Parker (1958). This model has 
the advantage of being fully controlled by only two parameters, 
namely the temperature Tp and the mass loss rate M, or equiv- 
alently Tp and the number density no at the base of the corona. 
The relation between these parameters writes 



M = 47tR 



2"(""P l^^^P l-^^^i/l-A 



OTp \ 4 



-re 



where R is the radius of the star and with 
GM/R 

IkTplm^ ' 

The velocity solution of the wind is implicitly given by 
„2,, A^R^ (3 AR] 



(30) 



(31) 



(32) 
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Study, Robinson et al. (2001) derive no = 1.7 x 10'^ cm ^ as an 
average for flare phases. We thus decided to assume as typical 
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Fig. 11. ySsw ratio as a function of the distance to the star, computed for 
a spherical silicate grain with radius .v = 0. 1 yum, in quiescent (thick 
solid line) and flare (thin solid line) states. The dotted line is the time- 
averaged ySsw, assuming that flares are present 10% of the time (y = 
0.1). 



where w - ^m^l{2kT^)Vsvi is the velocity normalized to the 
thermal velocity. With the above relations, we may compute 
the j6sw ratio as a function of the distance r from the star and 
of the parameters Tp and iiq. The analytical expression is rather 
complex and makes use of Lambert's W function that comes 
from the solution of Eq. (32), but we are interested here into 
the numerical values for/3sw. 

Our first remark is that contrary to 2pr, the coupling co- 
efficient Cd is here totally independent from the size of the 
grain. This means that for spherical grains, ySsw will be exactly 
oc i"' (see Eq. 28), for any grain size. This behaviour was noted 
by Mukai & Yamamoto (1982). Hence for sufficiently small 
grains, ySsw may reach arbitrarily high values. However, Mukai 
& Yamamoto (1982) show that in the case of the solar wind 
characterized by Tp ^ 10^ K and no - lO^cm"-', /3sw remains 
several orders of magnitude less than ySpr for usual grain sizes 
(ySsw would reach 0.5 or 1 for radii less than 1 nm). As a matter 
of fact, in the solar environment, the wind pressure is usually 
compared to Poynting-Robertson drag rather than to the full ra- 
diation pressure, with a ratio to PRD reaching ~ 0.3 (Gustafson 
1994). 

The situation may be drastically different with an active 
star like AUMic. Theoretical studies have shown that active 
M-dwarfs may have wind mass loss rates typically reaching a 
few times 10"'^ Mq yr"' (Wargelin & Drake 2001; Wood et al. 
2002), i.e. more than 100 times the solar value. Concerning 
AUMic explicitly, there are several density and temperature 
estimates in the coronal plasma near the star, deduced from 
the observation and analysis of spectral lines of highly ion- 
ized species. Monsignori-Fossi et al. (1996), based on the 
study of EUV lines of Fe xxi and Fe xxiv, give a temperature 
Tp ^ 8 X 10*" K in quiescent phases and 2-5 x 10^ K in flare 
phases; they also derive electron densities (that we will take 
equal to no) between 3 x lO'^ and 2 x lO'-' cm"-' during flares. 
Pagano et al. (2000) find no = 6.3x10'" cm"-' in quiescent state 
with some regions with no = 5 x 10" cm"^. Completing that 



values for quiescent state (no - 6.3x10 cm ,rp ^ 8x10 K), 
and (no = 1.7 x lO'^cm"^ Tp ^ 3 x 10^ K) during flares. This 
results via Eq. (30) in M = 9.4 X lO'^^Mgyr"' in quiescent 
state and a peak of M = 4.9 x 10"" Mq yr"' during flares. If 
we assume that flares are present 10% of the time, the average 
mass loss rate is 5.7 x 10"'^ Mq yr"' ^ 3 x lO^M©, which is in 
perfect agreement with the theoretical estimates quoted above. 
In particular. Wood et al. (2002) derives a power-law relation- 
ship between the X-ray activity of dwarfs and their mass loss 
rate. Applied to AUMic, this gives M ~ 1-2 x lO^'^Mgyr"' 
(Plavchan et al. 2005). 



6.3. Origin of the disk surface brightness profile 

We are now able to compute the j6sw ratio for any grain size 
at any distance from the star This is done in Fig. 11, where 
ySsw is plotted as a function of the stellar distance, in quiescent 
and flare regimes, for a grain of radius i = 0.1 /im. The cal- 
culation is done with Co computed from Eq. (27) neglecting 
sputtering (Cq turns out to remain between 2.0 and 2.1 over 
the whole range of stellar distances considered). The behaviour 
for other grain sizes may be directly deduced from this plot, as 
we strictly haveySsw '^ * ' . First of all, we note that there is only 
little variation of /3sw over a wide range of distances from the 
star. In fact, from Eq. (28) we see that the dependence in r of the 
ySsw ratio is entirely contained in the factor VswCd. Cd remains 
very close to 2, and Vsw is a moderately increasing function of 
r. ySsw is thus an increasing function of r, but the relative in- 
crease is only 12.4% between r = 15 AU and r = 100 AU, and 
8.4% between r = 100 AU and r = 500 AU. Hence considering 
ySsw as independent from r, like/3pr, is an acceptable approxima- 
tion. This justifies the introduction of the ySsw ratio, and shows 
that the dynamics of particles subject to wind pressure will be 
somewhat similar to that of particles subject to radiation pres- 
sure. In particular, we expect that grains will be pushed away 
in the outer disk as soon as ySsw will approach 0.5. 

Figure 1 1 shows that such large values of /3sw can easily be 
reached. For s - 0.1 /im, even in quiescent phase, ySsw - 0.4, 
the exact value depending on the assumed mass loss rate. But 
because of the recurrent flare episodes, we must consider the 
temporal average of y6sw. This is plotted in Fig. 11 as a dot- 
ted line, assuming that flares are present 10% of the time. Due 
to flares, grains with s = 0.1 /im are actually blown away by 
the stellar wind. As ySsw ^ i"', it turns out that all grains with 
s ^ lyum are able to diffuse into the outer disk, and even be 
ejected from the system. Actually the exact behaviour at large 
distance may depend on the location of the heliopause of the 
stellar wind, but this is not expected closer to a few hundreds 
of AU. 

The fact that ySsw is almost constant over several orders of 
magnitude of stellar distances, and thatySsw is able to reach suf- 
ficiently large values, show that particles subject to wind pres- 
sure will behave almost as if they were subject to a compara- 
ble radiation pressure. This explains why the surface brightness 
profile along the disk midplane (oc r"'* ' "^) resembles that of the 
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P Pictoris disk, the latter being due to the effects of the intense 
radiation pressure (Lecavelier des Etangs et al. 1996; Augereau 
et al. 2001; Thebault & Augereau 2005). 

From this study we conclude that in the AU Mic system, the 
stellar wind pressure is expected to play almost the same role as 
the radiation pressure in other debris disks around earlier type 
stars like yS Pictoris. Interestingly, Roberge et al. (2005) also 
came to the conclusion that the short lifetime of the circumstel- 
lar gas in the AU Mic system can only be explained by dissi- 
pation mechanisms other than photo-evaporation, including a 
possible stellar wind. Plavchan et al. (2005) also proposed that 
the dynamics of the grains in the inner AU Mic system could 
be controlled by the drag component of the stellar wind. Given 
that the collision time-scale we find within ~ 30 AU compares 
to the stellar wind drag time-scale inferred by Plavchan et al. 
(2005) (time-scale rescaled assuming M ^ 3 x lO^M©), the drag 
force may play a significant role in explaining the disk clearing 
in the inner regions. Interestingly, we find that the role of the 
flares, neglected in previous studies, is significant. 

7. Summary and Conclusion 

We have calculated the surface density profiles of the edge-on 
debris disk around the M-type star AU Mic by performing a 
direct inversion of two brightness profiles: the H-band profiles 
obtained by Liu (2004a) with the Keck adaptive optics system 
and the HST/ACS V-band profiles from Krist et al. (2005). We 
have solved the integral equation that links the scattered light 
surface brightness profile along the disk midplane to the disk 
surface density. The solution to this equation only depends on 
the light scattering phase function characterized by the asym- 
metry parameter \g\, which is constrained to be smaller than 
0.6. From the analysis of the H-band observations, we find that 
the dust distribution peaks around 35 AU, close to the location 
of the power law breakup in the surface brightness profile. The 
surface density profiles inferred from the inversion of the vis- 
ible images are quite different, and this difference may partly 
come from dust properties varying throughout the disk mid- 
plane. Even if the exact azimuthal profile of the disk cannot 
be unambiguously derived, the disk turns out to be obviously a 
non-symmetrical ring, as the two NE and SW branches are not 
identical. 

From the analysis of the disk colors measured in the visible 
with the HST/ACS instrument (Krist et al. 2005), we constrain 
the grain size distribution to be steeper than the classical s^^-^ 
differential size distribution. Grains as small as 0.1 /zm in ra- 
dius are in addition required to explain the observations when 
astronomical silicates are assumed. We infer a total mass of sil- 
icates grains smaller than 1 mm equal to 3.5-8 x 10""* M^. This 
dust mass compares the mass of the Ceres asteroid in the Solar 
System, but it is 14 to 30 times smaller than the dust mass in- 
ferred from millimeter observations by Liu et al. (2004b). It is 
also 9 to 20 times smaller than the mass we infer from the disk 
SED fitting. This difference cannot be explained by the disk 
optical thickness since the disk is optically thin in all direc- 
tions and at all wavelengths (xy < 0.04 at visible wavelengths). 
Possible explanations for this discrepancy include the fact that 
the simple power law size distribution that we find from the 



interpretation of the visible observations can probably not be 
extrapolated all the way to millimeter grain sizes (see e.g. the 
discussion in Thebault et al. 2003). It may also be possible that 
assuming hard spheres to calculate the grain optical properties 
is not appropriate if the grains have complex structures. 

We then explored whether the similarities between the sur- 
face brightness profiles of the AU Mic and fi Pictoris disks are 
coincidental or if they could result from similar dynamical pro- 
cesses. The global properties of the dust distribution in the 
[} Pictoris disk are very well explained by the role of radiation 
pressure. In this context, the breakup in surface brightness (ob- 
served at 120 AU) corresponds to the outer edge of the plan- 
etesimal disk that produces the dust particles (Augereau et al. 
2001). A very similar picture may apply to the AU Mic disk but 
at a distance of ~ 35 AU. The main difference lies in the dust 
diffusion mechanism. The radiation pressure from AU Mic ap- 
pears inefficient to diffuse dust into the outer disk even when 
the X-ray and UV flares are taken into account. Conversely, the 
corpuscular pressure from a stellar wind is very likely to be 
much more efficient. In fact, in the AU Mic environment, the 
wind pressure turns out to play an almost identical role as the 
radiation pressure in the fi Pictoris disk, with in both cases, a 
constant (or almost constant) ratio to stellar gravity. This ex- 
plains why the surface brightness profiles of the two disks are 
similar, despite different generating mechanisms. 

This comparison between debris disks about early and late- 
type stars would be generic if the additional role of flares was 
not significant in the case of AU Mic. The flares contribute to 
considerably enhance the effect of the stellar wind, and con- 
sequently to transport in the outer disk many additional parti- 
cles that would otherwise remain close to their parent bodies. 
Debris disks may be common around late-type stars, but an 
intense and frequent coronal activity may help to make them 
sufficiently extended to be detectable. 
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Appendix A: Numerical scheme to invert Eq. (5) 

The product integration method is a simple but efficient and 
generally numerically stable technique that can be used to in- 
vert an integral equation. This Appendix describes the main 
steps towards the implementation of this method in order to 
solve Eq. (5). We refer to Sect. 2 and 3.3 for the notations used 
in this appendix. 

We expand Eq. (5) as a sum of n integrals over ranges of 
equal size h 



T7^ Jr„,,-(i+l)h 



A(r,y) 

■ ylf2 _ y2 



dr 



where y - r^ax - nh and 

A{r, y) = O-D-' (r, 6) 0(r, 6) (f(0) + /(tt - 6)) 
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We note rj = r^^x - jh which gives r„ = y and r^ = r^^,^. S(y) 
and A(r, y) are supposed to be two known functions of r and/or 

y- 

We assume that for r within the range r^+i < r < rj, 
Cscti'^oir) A{r,y) can be approximated, to first order, by its 
value at the middle of the interval of size h (product integral 
method), namely at the position r - r^^i. The smaller is h, 
the better is the approximation. The surface brightness profile 
S{y) - S (r„) thus reads 



where 



and 



i=0 

I 



, n, r„) = I 

•Jrm 



•V^ 



■dr- arctan I 
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h being small compared to the size of the disk, taking n = 1 in 
the above expression of S{r„) leads to a first crude estimate of 
o"sca '^oif) close the outer edge r^^x of the disk 



0-sca2o(ri) 



Sin) 



A(ri,r„)I(ri,ro,r„) 
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For n > 2, we can write 
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S (r„) ^ 2_j '^sca 2o(r^+ 1 )A(r^._|_ 1 , r„)I(rj+u rj, r„) + 
o"sca2:o(r„_i )A(r„_i, r„)/(r„, r„_i, r„) 
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and then 

O'sca£o('-„_i) 



S(r„) 



A(r„_i,r„)/(r„,r„_i,r„) 



-.n-2 



2"=o O"sca2o(r^+i )A(r^^i , r„)/(rj+i, rj, r„) 
A(r„_i,r„)/(r„,r„_i,r„) 

" 2 

With this simple algorithm, we can directly invert Eq. (5) to get 
the product of the surface number density 'Loir) by the scatter- 
ing cross section osc-^ as a function of the distance to the star 
(Tfic-^ So is reconstructed from the outer disk edge (rmax) down to 
the closest possible distance. At each distance r from the star, 
o"sca So('") is thus calculated using the values of osc-i 2o previ- 
ously estimated at distances larger than r. We checked that with 
this simple algorithm we could reproduced the analytic results 
obtained by Nakano (1990) for an edge-on disk with a decreas- 
ing power-law surface brightness profile. 



Appendix B: Accounting for non-axisymmetric 
dislts: tlie @(r, B) smoothing function 

The inversion of the NW and SE surface brightness profiles do 
not depend on the position angle on the skyplane at which they 
has been measured. To simplify the discussion, and as far as 
the inversion is concerned, we will thus always assume that the 
brightness profile is measured along the >'-axis and orientated 
toward the positive y direction (see Fig. 1). 



If h (as defined in Appendix A) could be infinitely small, 
the inversion of Eq. (5) would lead to the reconstruction of the 
surface density of the disk precisely along the y-axis, namely 
at a polar angle 6 = jt/2. But because oscj^oif) is estimated 
at the middle of an interval of finite size h, the surface den- 
sity is calculated at polar angles 6o that slightly deviates from 
jt/2 and with lim Oq = 7r/2 when /z — > 0. The smoothing func- 
tion 0(r, ff) introduced in Sect. 3.3 in order to account for non- 
axisymmetric disks reads 



&SE(r,e) 



1 



00 j 2 



for the SE side of the disk, and 

e,,(,,).^x(-).^ 



(B.l) 



(B.2) 



for the NW side with < 6 < 7t/2 and where p is given by 

p = [o-sca2o(r)]Nw / [o-sca2o(r)]sE (B.3) 

Appendix C: Dynamics of particles with 
time-variable radiation pressure 

We investigate here the dynamics of particles with a time- 
variable ySp,. We assume that the star periodically alternates qui- 
escent and flare phases, respectively characterized for a given 
grain byySprq andySpr^ ratios, withySpr/ > ySpr.q- The star is as- 
sumed to spend a fixed fraction y of the time in flare phases. 
Strictly speaking, the grain follows bits of Keplerian orbits al- 
ternatively characterized by dynamical masses M{1 -/3pr,q) and 
M(l - y6pr,f). What is the global effect ? First, we should note 
that this problem cannot be treated by a classical perturbation 
method typical for celestial mechanics, as the "perturbation" 
(the flares effect) overcomes the main Keplerian forces. 

Let us consider a simple episode, characterized by a first 
quiescent phase (ySpr = /?pr,q), a flare (fip^ = y^pr.f), and a second 
quiescent phase (Jip^ = y6pr,q). During each of these phases the 
grain follows a different Keplerian orbit. Due to the flare, the 
Keplerian orbit during the first quiescent phase differs from that 
during the second quiescent phase. We want here to compute 
the orbital change. 

First, we note that during the whole episode, the angular 
momentum remains unchanged. Even if j6pr is time-variable, 
the resulting acceleration keeps pointing towards the central 
star, keeping the angular momentum constant. Hence the or- 
bital plane of the motion remains unchanged and so the quan- 
tity a(l - e^) during the successive quiescent phases. Second 
the orbital energy 6E change between the two quiescent phases 
may be computed easily, using the continuity of the velocity 
during the whole episode. After a simple algebra, we derive 



5E ^ GM(J3p,s - Pvm) 



nri 



(C.l) 



where ri and r2 are the distances of the particle to the star at 
the beginning of the flare and at the end respectively. Now, the 
important remark is that the period of recurrence of the flares 
(a few days at most) at far less than the orbital period of the 
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particle (hundreds of years at 35^0 AU). The distance change 
during the whole episode may be regarded as infinitesimal. 

Let us now consider the time span between the beginning 
of a flare and the beginning of the next flare. During this time 
span, we first have a flare that generates a (quiescent) orbital 
energy change described by Eq. (C.l), and a quiescent phase 
that keeps is constant up to the next flare. The time-span con- 
sidered may be regarded as infinitesimal compared to the or- 
bital period, so that the motion of the particle during this phase 
is nearly uniform. If dr is the infinitesimal distance change dur- 
ing this phase, we should have r2- r\ ^ ydr in Eq. (C.l). The 
orbital energy change within dr reads then 



d£ = GM(ySp,,f-/3p,,q)^. 



(C.2) 



Integrating this equation, we derive the macroscopic energy 
change that corresponds to the stellar distance change from r 
to r' : 



hE^GMy{p^,S-Pvm)\--\ 



This is equivalent to say that the quantity 

GMr(^p,,f-/3p,q) 

ieff - -fc H 

r 
Given the expression of the energy E, we derive 

1 2 GM(l-/?p,,q) GMy{l5^,s - Pvm) 
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1 2 GM(1 -;8p,,etf) 

= -V 

2 r 



(C.3) 



(C.4) 



where 



j6pr,eff = (1 - r)ySpi-,q + r/^pr.f • 



(C.5) 



(C.6) 



iieff is just the energy corresponding to a Keplerian orbit with 
central mass M(l -ySpr,etf), i-e. to a particle feeling a/Jpr = ySpr.eff 
ratio from the star Combined to the preservation of angular 
momentum, we see that on average, the dust particle will keep 
following a Keplerian orbit around the star, but characterized 
by ajSpr = ySpr.etf ratio. We checked this result numerically with 
a dedicated integrator. 

/?pr,efF is only the temporal average of the ySpr ratio between 
flare and quiescent phases. Finally, the net results is that when- 
ever the radiation pressure is time-variable, the grain behaves as 
if it was feeling a/Jpr ratio from the star that is just the temporal 
average of the variable y6pr. We should nevertheless not forget 
that this is valid if the recurrence frequency of the flares largely 
overcomes the orbital frequency of the orbit. Significantly dif- 
ferent outcomes should be expected for instance if there was a 
resonance between the mean motion of the orbit and the recur- 
rence frequency of flares, which is not the case here. 
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